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We report the observation of B^-B° oscillations from a time-dependent measurement of the B®-B® 
oscillation frequency Am,. Using a data sample of 1 fb" 1 of pp collisions at y/s = 1.96 TeV col- 
lected with the CDF II detector at the Fermilab Tevatron, we find signals of 5600 fully reconstructed 
hadronic B s decays, 3100 partially reconstructed hadronic B s decays, and 61 500 partially recon- 
structed semileptonic B a decays. We measure the probability as a function of proper decay time that 
the B a decays with the same, or opposite, flavor as the flavor at production, and we find a signal for 
Bg-Bs oscillations. The probability that random fluctuations could produce a comparable signal is 
8 x 10~ 8 , which exceeds 5<r significance. We measure Am s ~ 17.77 ±0.10 (stat) ±0.07 (syst)ps -1 
and extract \V t d/V ts \ = 0.2060 ±0.0007 (exp) to'.ooeo (theor). 



PACS numbers: 14.40.Nd, 12.15.Ff, 12.15.Hh, 13. 20. He 



Since the first observation of particle-antiparticle transformations in neutral B mesons in 1987 |l|], the determination 
of the Bg-Bg oscillation frequency Am s from a time-dependent measurement of B®-B® oscillations has been a major 
objective of experimental particle physics 0]. This frequency can be used to extract the magnitude of Vt a , one of the 
nine elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix Q . Recently, we reported Q the strongest evidence 
to date of the direct observation of B®-B® oscillations. That analysis used 1 fb -1 of data collected with the CDF II 
detector at the Fermilab Tevatron, and the probability that random fluctuations would produce a comparable signal 
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was 0.2%, corresponding to 3<r signal significance. This level of significance is insufficient to claim a firm observation, 
however, under the oscillation hypothesis we determined Am s — 17.31 tp'is ( s t & t) ±0.07 (syst)ps -1 . In this Letter 
we report an update of this measurement that uses the same data set with an improved analysis and reduces this 
probability to 8 x 10~ 8 (> 5c), yielding the first definitive observation of time-dependent B®-B® oscillations. 

We improve the analysis in Ref. |4| by increasing the B s signal yield and improving the performance of the methods 
used to identify the flavor (b or b) of the B s at production. The previous analysis used B s decays in hadronic 
(S° — » _D+7r~, _D+7r~7r + 7r~) and semileptonic (B° — > Df^l^Di, I = e or /x) decay modes 0- We used Df — ► 4>ir + , 
K*(892) Q K + , and ir + TT~ir + , with (j> — > K + K~ and A'* — » X~7r + . Several improvements lead to increased signal 
yields. We use particle identification techniques to find kaons from D s meson decays, allowing us to relax kinematic 
selection requirements on the D s decay products. This results in increased efficiency for reconstructing the D s 
while maintaining excellent signal to background. In the hadronic channels, we employ an artificial neural network 
(ANN) to improve candidate selection resulting in larger signal yields at similar or smaller background levels. The 
ANN selection makes it possible to use the additional decay sequence B® — ^ D+ix~it + it~ , with Df — > 7r + 7r~7r + , as 
well. We add significant statistics using partially reconstructed hadronic decays in which a photon or ir° is missing: 
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Finally ANNs are used to enhance 



the performance of the methods used to identify the flavor of the B s at production. With all these improvements, the 
effective statistical size of our data sample is increased by a factor of 2.5. 

The CDF II detector [|| consists of a magnetic spectrometer surrounded by electromagnetic and hadronic calorime- 
ters and muon detectors. The key components for this measurement are listed below. Precision determination of the 
decay point is provided by a seven-layer double-sided silicon-strip detector and a single-sided layer of silicon mounted 
directly on the beampipe at an average radius of 1.5 cm. A 96-layer drift chamber is used for both precision tracking 
and dE/dx particle identification. Time-of- flight (TOF) counters surrounding the drift chamber are used to identify 
low-momentum charged kaons. A three-level trigger system selects, in real time, pp collisions containing charm and 
bottom hadrons by exploiting the kinematics of production and decay, and the long lifetimes of D and B mesons. A 
crucial component of the trigger system for this measurement is the Silicon Vertex Trigger, which makes it possible 
to collect our large sample of B s mesons in the fully or partially reconstructed hadronic decay modes, giving CDF 
unique sensitivity to B s oscillations. 
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FIG. 1: (Left panel) The invariant mass distributions for the D~^{(fm + ) candidates [inset] and the £~ Df (tf>ir + ) pairs. The 
contribution labelled "false lepton & physics" refers to backgrounds from hadrons mimicking the lepton signature combined 
with real D s mesons and physics backgrounds such as B° — > DfD~, Dt — > 4>tt + ,D~ £~X. (Right panel) The invariant 
mass distribution for 3° —> Df (cj>Tv + )iY~ decays including the contributions from — » D* s + -k~ and — > Dt p~ ■ In this 
panel, signal contributions are drawn added on top of the combinatorial background. 
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To reconstruct B® candidates, we first select D+ candidates. These Df candidates are combined with one or three 
additional charged particles to form Df£~, D+7r~, or D+7r~7r + 7r~ candidates. In the previous analysis, we reduced 
combinatorial backgrounds by applying requirements on selection quantities such as the minimum pj- |7j of the B® and 
its decay products, and the quality of the reconstructed B® and Dj decay points and their displacement from the pp 
collision position. In this analysis, we add a kaon identification likelihood formed from TOF and dE/dx information. 
For decay modes with kaons in the final state, we use this likelihood to reduce combinatorial background from random 
pions or physics backgrounds such as D + — » A"~7r + 7r + . In JJ, we vetoed Df candidates consistent with the D + mass 
hypothesis, which resulted in a substantial loss of signal efficiency. Kaon identification makes it possible to relax 
kinematic requirements (charged particle pt and the D + veto) leading to a substantial increase in signal efficiency. 

In the semileptonic channel, the main gain is in the Df£~, Df — > K*(892)°K + sequence, where the signal is 
increased by a factor of 2.2. An additional gain in signal by a factor of 1.3 with respect to our previous analysis comes 
from adding data selected with different trigger requirements. In total the signal of 37000 semileptonic B s decays in 
is increased to 61 500, and the signal to background improves by a factor of two in the sequences with kaons in 
the final state. The distributions of the invariant masses of the Df(<pTT + )£~ pairs mr) s t and the Df(<pTr + ) candidates 



are shown in Fig. ^ We use mD 3 i to help distinguish signal, which occurs at higher mu s 
physics (e.g., double-charm decays of B mesons) backgrounds. 
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data (a) 
fit 

B^D;(i?K)7i 
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□ B°-> D + 7T 
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— fit 
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FIG. 2: The invariant mass distributions for B° — > D^n~ (top panels) and Dfir~Tr + -K~ (bottom panels). Signal contributions 
are added on top of the combinatorial background. Contributions from partially reconstructed B a decays are taken into account 
in the fit and are not shown. 



In the hadronic decay modes, we use an ANN to enhance the signal selection of the previous analysis. The ANN 
uses quantities such as the selection criteria listed above as well as the kaon identification likelihood. The network is 
trained using simulated signals generated with Monte Carlo methods. For combinatorial background, we use sideband 
regions in the mass distribution of the B s candidates from data. In this analysis, we add the partially reconstructed 
signal between 5.0 and 5.3GeV/c 2 from B° a — > Z?*+7r~, D* + — > Dfj/n in which a photon or 7r° from the D* + is 
missing and B® — > Df p~ , p~ — > 7r~7r° in which a tt° is missing. The mass distributions for B® — > Dfir~, Df — > (f)ir + 
and the partially reconstructed signals are shown in Fig. ^ The mass distributions for the other five hadronic decay 
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Decay Sequence Signal S/B gain w.r.t \4j 

TO 13% 

3.4 n.a. 

2.0 35% 

2.1 22% 
2.7 92% 
1.1 110% 
2.6 n.a. 



TABLE I: Signal yields (S) and signal to background ratio (S/B) in the various hadronic decay sequences. The gain refers to 
the percentage increase in S/sJS + B relative to Q. 



sequences are shown in Fig. [21 In these modes, we require the masses of the candidates to be greater than 5.3 GeV/ c 2 . 
Candidates with masses greater than 5.5GeV/c 2 are used to construct probability density functions (PDFs) for 
combinatorial background. Table [I] summarizes the signal yields. 

The reconstructed decay time in the B s rest frame is t = tob^t/Pt 00 * 1 ' where Lt is the displacement of the B s 
decay point with respect to the primary vertex projected onto the B s transverse momentum vector, and p™ con is the 
transverse momentum of the reconstructed decay products. In the semileptonic and partially reconstructed hadronic 
decays, we correct t by a factor k = p r rp con /pt(B s ) determined with Monte Carlo simulation (Fig. EJl ■ 

The decay time resolution 174 has contributions from the momentum of missing decay products (due to the spread 
of the distribution of k) and from the uncertainty on Lt- The uncertainty due to the missing momentum increases 
with proper decay time and is an important contribution to a t in the semileptonic decays. To reduce this contribution 
and make optimal use of the semileptonic decays, we determine the n distribution as a function of Tnu s i (Fig. EJ. We 
estimate the contribution from the uncertainty on Lt to a t for each candidate using the measured track parameters 
and their estimated uncertainties. 

The distribution of at for fully reconstructed decays has an average value of 87 fs, which corresponds to one fourth 
of an oscillation period at Am s — 17.8 ps _1 . The distribution is nearly Gaussian with an rms width of 31 fs. For the 
partially reconstructed hadronic decays, the average at is 97 fs, and the addition to at due to the missing photon or 
ir° is very small (Fig. |2J. For semileptonic decays, at is worse due to decay topology and the much larger missing 
momentum of decay products that were not reconstructed. The increase of at with t is illustrated in Fig.|3for different 
ranges of mu 3 i- 

The flavor of the B s at production is determined using both opposite-side and same-side flavor tagging techniques. 
The effectiveness Q = eT> 2 of these techniques is quantified with an efficiency e, the fraction of signal candidates with 
a flavor tag, and a dilution V = 1 — 2w, where w is the probability that the tag is incorrect. 

At the Tevatron, the dominant 6-quark production mechanisms produce bb pairs. Opposite-side tags infer the 
production flavor of the B s from the decay products of the b hadron produced from the other b quark in the event. 
In the previous analysis, we used lepton (e and /1) charge and jet charge as tags, and if both types of tag were 
present, we used the lepton tag. In this improved analysis, we add an opposite-side flavor tag based on the charge 
of identified kaons, and we combine the information from the kaon, lepton, and jet charge tags using an ANN. The 
dilution is measured in data |8| using large samples of B~ , which do not change flavor, and B°, which can be used 
after accounting for their well-known oscillation frequency. The combined opposite-side tag effectiveness improves by 
20% to Q — 1.8 ± 0.1 %. Most of the improvement is for candidates with both a lepton and jet-charge tag. 

Same-side flavor tags are based on the charges of associated particles produced in the fragmentation of the b quark 
that produces the reconstructed B s . In the previous analysis, we used a same-side tag based on our kaon particle- 
identification likelihood; here we use an ANN to combine our kaon particle-identification likelihood with kinematic 
quantities of the kaon candidate into a single tagging variable T. Tracks close in phase space to the B s candidate 
are considered as same-side kaon tag candidates, and the track with the largest value of T is selected as the tagging 
track. We predict the dilution of the same-side tag using simulated data samples generated with the PYTHIA Monte 
Carlo program. The predicted fractional gain in Q from using the ANN is 10%. Control samples of B~ and B° 
are used to validate the predictions of the simulation. The effectiveness of this flavor tag increases with the pr of 
the B®\ we find Q = 3.7% (4.8%) in the hadronic (semileptonic) decay sample. The fractional uncertainty on Q is 
approximately 25% Q. If both a same-side tag and an opposite-side tag are present, we combine the information 
from both tags assuming they are independent. 

We use an unbinned maximum likelihood fit to search for B s oscillations. The likelihood combines mass, decay 



Bs ' -> DJ{<piT t )tt- 2000 

Partially reconstructed 3100 

Bs -> D+(K* (892)° K + )tt- 1400 

Bs -> Df (TT + lT-TT + )lT- 700 

B° -> D+(^7r+)7r"7r + 7r- 700 
13° ~* D+(K* (892)° K + )ty--r + ty~ 600 

B° — > D+(7r + 7r~7r + W~7r + 7r~ 200 
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FIG. 3: (Left panel) The distribution of the correction factor n in semileptonic and partially reconstructed hadronic decays 
from Monte Carlo simulation. (Right panel) The average proper decay time resolution for B s decays as a function of proper 
decay time. 



time, decay-time resolution, and flavor tagging information for each candidate, and includes terms for signal and each 
type of background. Details of the fit are described in [iHloj]. 

Following the method described in ^l|, we fit for the oscillation amplitude A while fixing Am s to a probe value. The 
oscillation amplitude is expected to be consistent with A = 1 when the probe value is the true oscillation frequency, 
and consistent with .4 = when the probe value is far from the true oscillation frequency. Figure 01 shows the fitted 
value of the amplitude as a function of the oscillation frequency for the semileptonic candidates alone, the hadronic 
candidates alone, and the combination. The sensitivity |j,[ll( is 19.3ps _1 for the semileptonic decays alone, 30.7PS" 1 
for the hadronic decays alone, and 31.3 ps -1 for all decays combined. At Am s — 17.75 ps -1 , the observed amplitude 
A = 1.21 ± 0.20 (stat.) is consistent with unity, indicating that the data are compatible with B S -B S oscillations with 
that frequency, while the amplitude is inconsistent with zero: Aj oa — 6.05, where is the statistical uncertainty 
on A (the ratio has negligible systematic uncertainties). The small uncertainty on A at Am s = 17.75 ps -1 is due to 
the superior decay-time resolution of the hadronic decay modes. 

We evaluate the significance of the signal using A = log{£ A=0 / £^ =1 (Am s )], which is the logarithm of the ratio 
of likelihoods for the hypothesis of oscillations (A — 1) at the probe value and the hypothesis that .4 = 0, which is 
equivalent to random production flavor tags. Figure 01 shows A as a function of Am s . Separate curves are shown for 
the semileptonic data alone (dashed), the hadronic data alone (light solid), and the combined data (dark solid). At 
the minimum Am s = 17.77 ps^ 1 , A = —17.26. The significance of the signal is the probability that randomly tagged 
data would produce a value of A lower than —17.26 at any value of Am s . We repeat the likelihood scan 350 million 
times with random tagging decisions; 28 of these scans have A < —17.26, corresponding to a probability of 8 x 10~ 8 
(5.4 a), well below 5.7 x 10" 7 (5 a). 

To measure Am s , we fix A — 1 and fit for the oscillation frequency. We find Am s = 17.77 ± 0.10 (stat) ± 
0.07 (syst) ps" 1 . The only non-negligible systematic uncertainty on Am s is from the uncertainty on the absolute scale 
of the decay-time measurement. Contributions to this uncertainty include biases in the primary- vertex reconstruction 
due to the presence of the opposite-side b hadron, uncertainties in the silicon-detector alignment, and biases in track 
fitting. The uncertainty on the correction k for the hadronic candidates with a missing photon or 7r° is included and 
has a negligible effect. 

The B®-B® oscillations are depicted in Fig. [5j Candidates in the hadronic sample are collected in five bins of proper 
decay time modulo the measured oscillation period 2n/Am s . In each bin, we fit for an amplitude (the points in 
Fig. EJ) using the likelihood function .4j, which takes into account the effects of background, flavor tag dilution and 
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FIG. 4: The measured amplitude values and uncertainties versus the B®-B® oscillation frequency Am s . (Upper Left) Semilep- 
tonic decays only. (Lower Left) Hadronic decays only. (Upper Right) All decay modes combined. (Lower Right) The logarithm 
of the ratio of likelihoods for amplitude equal to one and amplitude equal to zero, A = log[£ A=0 / £. A=1 (Am,)], versus the 
oscillation frequency. The horizontal line indicates the value A = —15 that corresponds to a probability of 5.7 x 10~ 7 (5<r) in 
the case of randomly tagged data. 



decay-time resolution for each candidate. The curve shown in Fig. [5] is a cosine with an amplitude of 1.28, which is 
the observed value in the amplitude scan for the hadronic sample at Am s — 17.77ps _1 . As expected, the data are 
well represented by the curve. 

The measured B®-B® oscillation frequency is used to derive the ratio |Vfd/Vt s | = £y Sm^m"^ E3' ^ s m P u ^ s we 
use m B o/m B v = 0.98390 [13 with negligible uncertainty, Am d = 0.507 ± 0.005 ps" 1 [l| and £ = 1.21 to.^l We 
find \V td /V ts \= 0.2060 ±0.0007 (exp) (theor). _ 

In conclusion, we report the first observation of B®-B® oscillations from a decay-time dependent measurement of 
Am s . Our signal exceeds 5<r significance and yields a precise value of Am s , which is consistent with standard model 
expectations. This result supersedes our previous measurement 0. 
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